On the nature of tunable hole g-factors in quantum dots 
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Electrically tunable g-factors in quantum dots are highly desirable for applications in quantum 
computing and spintronics. We report giant modulation of the hole g-factor in a SiGe nanocrystal 
when an electric field is applied to the nanocrystal along its growth direction. We derive a contribu- 
tion to the g-factor that stems from an orbital effect of the magnetic field, which lifts the Kramers 
degeneracy in the nanocrystal by altering the mixing between the heavy and the light holes. We 
show that the relative displacement between the heavy- and light-hole wave functions, occurring 
upon application of the electric field, has an effect on the mixing strength and leads to a strong 
non-monotonic modulation of the g-factor. Despite intensive studies of the g-factor since the late 
50's, this mechanism of g-factor control has been largely overlooked in the literature. 



In the past decade, a great effort has been devoted 
to the realization of spin qubits in semiconductors^. 
Spin manipulation was achieved through different ap- 
proaches: magnctic-ficld-driven electron spin resonance 3 -, 
electric-dipole spin resonance^—, and fast control of the 
exchange coupling^. Another possibility for electric- 
field spin manipulation is the g-tensor modulation reso- 
nance, which has been used on ensembles of spins in two- 
dimensional (2D) electron systems^. This technique re- 
lies on anisotropic and electrically tunable g-factors. Re- 
cently, several experiments have addressed the g-factor 
modulation by means of external electric field o 10 ' 11 , and 
different mechanisms were evoked to explain the observed 
g-factor tunability, such as compositional gradients^ and 
quenching of the angular moment u m 1 1 1 1 2 . Here we report 
the experimental observation of an exceptionally large 
and non-monotonic electric-field modulation of the hole 
g-factor in SiGe QDs. To interpret this finding we have 
to invoke a new mechanism that applies to hole-type low- 
dimensional systems. This mechanism relies on the ex- 
istence of an important, yet overlooked correction term 
in the g-factor whose magnitude depends on the mixing 
of heavy and light holes. We show that in SiGe self- 
assembled QDs an electric field applied along the growth 
axis can be used to efficiently alter this mixing and pro- 
duce large variations in the hole g-factor. 

Our SiGe QDs were grown by molecular-beam epi- 
taxy on a silicon-on-insulator substrate. The Stranski- 
Krastanow growth mode was tuned to yield dome-shaped 
QDs with height w = 20 nm and base diameter d = 
80 nm. A sketch of the device is shown in Fig. Q] (a) . 
The QD is contacted by two 20-nm-thick Al electrodes, 
acting as source and drain leads. A Cr/Au gate electrode 
is fabricated on top of the QD with a 6-nm-thick hafnia 
interlayer deposited by atomic-layer deposition. This top 
gate, together with the degenerately-doped Si back gate, 
allows a perpendicular electric field to be applied while 
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FIG. 1: (a) Schematic cross section of a SiGe QD device, (b)- 
(d) Color plot of dI Bd /dV sd (V tg , V sd ) for B z = 70 mT, 3T and 
5T, respectively (Vb g = 0). The lines indicated by rhombis 
correspond to the onset of tunneling via Zeeman-split levels 
for N — 1 and N + l holes on the QD. The lines indicated by 
a star and by a circle correspond to singlet-triplet excitations 
for N holes. 



maintaining a constant number of holes in the SiGe QD. 

Measurements of the g-factor were performed using 
single-hole tunneling spectroscopy. A typical differential 
conductance (dI sc i/dV s d) measurement as a function of 
top-gate voltage (Vt g ) and source-drain bias voltage (V s d) 
is shown in Fig. [T] (b). All measurements reported here 
were done in a 3 He refrigerator with a base temperature 
of 250 mK. In order to suppress the superconductivity 
of the leads, a small magnetic field, B z = 70 mT, was 
applied along the z axis, i.e. perpendicular to the (x, y) 
growth plane. Diamond-shaped regions, where the cur- 
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rent vanishes due to Coulomb blockade, can be clearly 
observed in Fig. [T] (b). The charging energy is about 
lOmeV. Outside the diamonds, additional lines denoting 
transport through excited orbital states can be observed. 
Figs. [1] (c) and (d) show the same Coulomb-blockade 
regime for B z = 3T and B z = 5T, respectively. The 
magnetic field causes a splitting of the diamond edges as 
indicated by green rhombis. This splitting follows from 
the lifting of Kramers degeneracy in the ground states 
associated with the side diamonds. We thus conclude 
that the central diamond corresponds to an even num- 
ber, N, of confined holesi. The Zeeman energy splitting 
is given by Ez = g±fiBB z , where [Ib is the Bohr magne- 
ton and g± is the absolute value of the g-factor along z. 
From the splitting of iV-hole diamond edges we extract 
g± = (3.0±0.4) &ndg ± = (2.8±0.4) for the TV- 1 and the 
N + 1 ground states, respectively. The line indicated by 
a star in Fig. Q] (b) is due to the spin-triplet excited state 
for N holes on the QD. We measure a 2meV singlet- 
triplet energy in this particular QD, which is an order 
of magnitude larger than for electrons in Si/SiGe het- 
erostructuresii We note that large singlet-triplet excita- 
tion energies are particularly desirable for the observation 
of spin blockade in double-dot experiments^. Upon in- 
creasing B z , the line denoted by a star splits as shown by 
the emergence of second parallel line, denoted by a circle, 
that shifts away proportionally to B z (see Figs.[T](c) and 
(d)). This behavior corresponds to the Zeeman splitting 
of the excited spin-triplet stated with g± = (2.8 ± 0.4). 
Hereafter, we will concentrate on g-factor measurements 
in spin-1/2 ground states. 

Our dual-gate devices allow us to measure the depen- 
dence of the g-factor on a perpendicular electric field, 
F, at constant number of holes. The principle of such 
a measurement is illustrated in Fig. [5] (a). The Zeeman 
splitting is given by the distance between the blue and 
the red circles along Vb g , multiplied by a calibration fac- 
tor a. The latter is obtained by dividing V s d by the 
distance between the green and the red circles. In order 
to investigate the F-dependence of the g-factor, a con- 
stant V s d = 2.6 mV was applied and Vb g was swept while 
stepping Vt g . The magnetic field was fixed at 4T. The 
data is shown in Fig. [2] (b) and the extracted g-factors 
are displayed in Fig. [U (c). We observe an exceptionally 
large g-factor modulation (Sg/g ~ 1) denoting a strong 
effect of the applied F. The g-factor increases slowly to 
a maximum value of 2.6 and then drops rapidly down to 
a point where the Zeeman splitting can no longer be re- 
solved. Comparably large g-factor variations have been 
observed in other similar measurements, see Appendix El 

In order to uncover the origin of this unusual behav- 
ior, we modelled the QD electronic states in terms of 
heavy-hole (HH) and light-hole (LH) subbands. Given 
the relatively large anisotropy of dome-shaped QDs, we 
initially considered the two-dimensional (2D) limit re- 
sulting from confinement along the growth axis. Due to 
quantum confinement along z = [001] and strain, the 4- 
fold degeneracy of the valence band at T-point is lifted. 
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FIG. 2: (a) Left: Color plots of dI sd /dV sd {V bg ,V sd ) for 
B - 70 mT and 4T. At 4T the Zeeman splitting is clearly 
visible. Right: Corresponding schematic diagram illustrat- 
ing the measurement principle to extract the Zeeman energy 
splitting (and hence the g-factor) from gate-voltage sweeps at 
constant V s d (see the horizontal green line) . (b) Color plots of 
dI a d/dV s d(Vbg, Vtg) for a fixed V s d = 2.6mV. These data sets 
demonstrate the modulation of g± by a perpendicular electric 
field proportional to Vt g — Vt g . (c) g±(Vb g , Vt g ) as extracted 
from (b). Below g± ~ 0.75 the Zeeman splitting cannot be re- 
solved any more due to the finite broadening of the tunneling 
resonances. 

The top-most subband has HH character and its in-plane 
dispersion relation is described by the effective 2D Hamil- 
tonian 

HcS = ( k l + + \9\\VB ((T X B X + (TyBy) 

-\^g^BO z B z + U{x,y), (1) 

where k x and k y are the in-plane momentum operators, 
m|| = m/ (71+72) is the in-plane effective mass^, gu = 3q 
and g± = 6k + ^-q are, respectively, the in-plane and 
transverse g-factor a 15 ' 16 , <x are the Pauli matrices in the 
pseudospin spaced, and U(x, y) is the in-plane confining 
potential in the QD. We use standard notations for the 
Luttinger parameters 71, 72, 73, k, and q^. Since g«K, 
it is appropriate to assume g± ~ 6k. The minus sign in 
front of \g_\_ in Eq. ([!} is introduced for the convenience 
of having g± positive for Ge. 

First we consider the possibility that the observed g- 
factor modulation arises from a compositional gradient. 

It is well known that Si and Ge intermix leading to 
the formation of a Sii_ x Ge x alloy in which x increases 
motonically with 2, being zero at the base (z = —w) and 
approaching unity at the apex (z = 0) of the QDi£. Since 
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Ksi = —0.42 and kq — 3.41. one would expect that g± 
increases with F following a vertical shift of the HH wave 
function towards the apex. This compositional- gradient 
mechanism was exploited in Al x Gai_ x As quantum wells 
to implement electrical control of electron spins&£. While 
for electrons this may well be the only efficient way to 
control the g-factor, the situation might be different for 
holes. 

To find an upper bound for the g-factor variation re- 
sulting from the compositional gradient, we take the 
steepest dependence reported for the Ge content across 
the QD±2, 



x(z) 



^max \ / ± i ; 

W 



-w < z < 0, 



(2) 



where w = 20 nm is the height of the dome-shaped QD. 
To account for the existing uniaxial strain, we assume 
that the in-plane lattice constant on increases linearly 
from 5.47 A at the base to 5.59 A at the apexi^. With 
these two ingredients, the valence band profiles E v (z) 
for all types of holes are calculated using interpolation 
schemes devised for SiG o 20 ' 21 (see inset of Fig. [3]). The 
HH ground state is thus confined to a triangular potential 
well arising from the compositional gradient. An electric 
field applied along z adds a term —eFz to E v (z). For 
a given F, the HH wave function ip{z) is obtained by 
solving the Schrodingcr equation numerically. The HH 
g-factor is found as a weighted average 



g± w 6 (k) = 6 / k [x(z)j \ip(z)\ 2 dz, 



(3) 



where k (x) is obtained using the non-linear interpolation 
described in Ref. [22|. The resulting dependence g±(F) 
is shown in Fig. [3] We distinguish two regimes: that 
of a strongly asymmetric (triangular) potential well and 
that of a symmetric potential well. The modulation of 
the g-factor is largest in the latter regime (see dotted 
line in Fig. |3J}; we obtain dg±/dF « 0.41 m/MV. While 
the magnitude of the modulation is close to what is ob- 
served in the experiment, the sign of dg±/dF is never- 
theless opposite (see Fig. [21 (c)). We conclude that the 
compositional-gradient cannot explain our data. There- 
fore, from now on, we shall discard this mechanism and 
assume the Ge content to be constant within the QD. 

We revisit the derivation of Eq. (JXJ) , starting from 
the Luttinger Hamiltonian. In the 2D limit, the 4x4 
Luttingcr Hamiltonian separates into 2x2 HH and LH 
blocks, see Appendix [Bl To leading order of w/d <C 1, 
the HH and LH sectors are connected by the off-diagonal 
mixing blocks 2 ^ 

H h i = (Hih)^ = » ^ 73 {k x a y + k y a x ) k z , (4) 

171 

where k x and k y are 2D versions of momentum operators 
(insensitive to in-plane magnetic fields), k z = —ihd/dz, 
and a x and o~ y are the Pauli matrices in a pseudospin 
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FIG. 3: Expected electric-field dependence of g± for a SiGe 
QD with a strong compositional gradient. The numerical re- 
sult (solid line) has two regimes, designated by fits to two 
simplistic models. The dashed line shows a fit to the expres- 
sion (ft) = Koo — Aft (1 + -F/Fintr) -1 '' 3 , derived for a trian- 
gular potential well with an intrinsic electric field -Fintr for 
z < and infinite barrier for z > 0. The dotted line shows 
a fit to a linear dependence, obtained for a symmetric po- 
tential well. In the latter regime, g± is most sensitive to F, 
with dg±/dF — 0.41 m/MV. At large negative F, the wave 
function is pushed into the Si-rich region, where g± becomes 
negative. (Inset) Energy profiles for the heavy-hole (HH), 
light-hole (LH), and split-off (SO) bands for a Ge content x 
given by Eq. ([2| with x max = 0.8. We complemented the QD 
model with a thin layer of Si substrate at z G [—24, —20] nm 
and a strained Si capping layer at z G [0, 4] nm. Strain is 
taken into account resulting in a splitting between HH and 
LH bands. 



space, introduced simultaneously for heavy and light 
holes 2 ^. 

The mixing blocks in Eq. (U) are proportional to k z . In 
spite of the fact that k z averages to zero for each type of 
hole separately, it cannot be discarded in Eq. because 
matrix elements of the type (tph \ k z \tpi) are, in general, 
non-zero and scale as 1/w for w — > 0. Here, i/jh(z) and 
ipi{z) obey two separate Schrodingcr equations, for heavy 
and light holes, respectively (see below). This observa- 
tion allows us to anticipate that in second-order pertur- 
bation theory the mixing blocks lead to an energy cor- 
rection containing H^Hih cx k z in the numerator and 
Hu — Hhh oc fc 2 in the denominator. This correction 
does not vanish in the 2D limit (k z — > oo). At the 
same time, the correction to the wave function vanishes 
as k\\/k z ~ w/d. 

Using second-order perturbation theory, we recover 
Eq. |T]) for the top-most hole subband. Yet, at the lead- 
ing (zeroth) order of w/d <C 1, we obtain the following 
modified expressions for the effective mass and perpen- 
dicular g-factor, 

27 



m\\ = : , 9± =6k+— q-2~f h . (5) 

7i+72-7/i 2 

The in-plane g-factor remains unchanged (gu — 3q) at 
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this order. In Eq. 7/, is a dimensionlcss parameter 
sensitive to the form of the confinement along z, 



M 

TO 



E 



[Ml*iM>l 

EL - E' 1 



(6) 



Here, the sum runs over the LH subbands and the wave 
functions tpn^ l (z) and energies E^ 1 obey 



2m 



where m^' = 771/(71 =F 272) and V/j/;(z) is the confining 
potential seen by the heavy/light hole. 

When Vh{z) and Vi(z) are infinite square wells, the 
quantity 7^ in Eq. (j6]) can be derived analytically, 



h/l 



+ V h/l (z) 



^J l (z)^E^/ l (z), (7) 



lh 



1271 
71 + 2 72 



1-/3 7T (1 - /?) 



(8) 



where /3 = m l A _/m']_+5E m i/E[, with <5-E oi = Vh — Vt be- 
ing the splitting of the valence band due to uniaxial strain 
and E[ = n 2 h 2 /2m l ± w 2 . Notably, one has ipn( z ) = *Pli( z ) 
in this case, because the masses and m l ± drop out of 
the expressions for the wave functions. With the appli- 
cation of an electric field, the wave functions 4>n an d V'i 
begin to shift relative to each other, because of their dif- 
ferent effective masses. Although jh can only be numer- 
ically computed, its qualitative dependence on F can be 
inferred by inspecting Eq. ([5]). Note that the 71 = 1 term 
in the sum has the smallest energy denominator and, 
therefore, it is expected to have a dominant contribution. 
For a square-well potential, however, this term vanishes 
by symmetry. As a result, the symmetric point F = 
corresponds to a minimum in 7/i(F), since E n > E\. 
Away from F = 0, jh increases quadratically, 7^ cx F 2 , 
up to the point where the electric field is strong enough 
to shift the HH wave function (eFw ~ E\ - E\). Then, 
7/1 increases roughly linearly up to the point where the 
LH wave functions begin to shift (eFw ~ E l 2 — E[). Upon 
further increasing F, 7^ increases weakly and saturates 
to a constant. We remark that g±_ is modified by 7^ even 
at feii = 0, despite the fact that no HH-LH mixing occurs 
at feii = 0. In fact, g^ is sensitive to orbital motion in 
a perpendicular magnetic field 2 ^, and even a small B z 
translates to fc|| ^ 0, leading to HH-LH mixing. 

Our result in Eq. represents the zeroth-order term 
in the expansion g = g^ +g^ + . . ., where g^ cx (w/d) 2 
is the subleading-order term. Unlike the main term, the 
correction g^ 2 > is sensitive to the in-plane confining po- 
tential U(x,y) and it originates from the HH-LH inter- 
ference terms in the wave function. We shall address g^ 2 ' 
in a separate work. In Fig. 21 we fit the experimental 
data using only the leading, zero-order term. The insets 
in Fig. @] illustrate how the HH and LH wave functions, 
ipi(z) (red) and ipx(z) (blue), shift upon application of 
the electric field. Note that the transition from the square 
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FIG. 4: The g±(F) dependence according to Eq. (JSJ (red line) 
and corresponding experimental data (open dots). (Insets) 
Schematics of the HH and LH confinement at different F. 
The central inset shows the symmetric-well configuration at 
F = 0. The upper insets show the regime of strong F where 
both types of holes are confined to one side of the (triangular) 
potential well. The lower insets show the intermediate regime 
in which F is strong enough to push the heavy hole to one 
side of the well, while the light hole remains spread over the 
entire well width. 



well (central inset) to the triangular well (upper insets) 
occurs in two steps. First, the HH wave function shifts by 
8z ~ w, while the LH wave function remains nearly unaf- 
fected (lowest insets) . Then, the LH wave function shifts 
as well (highest insets). At even larger F (not shown) g± 
saturates to g± ~ 0.6. The evolution of g± taking into 
account jh , shown as a solid red line in Fig. |4l qualita- 
tively reproduces our experimental results (open dots in 
Fig. 2]). 

Finally, we remark that the correct 2D limit of the Lut- 
tingcr Hamiltonian has been largely overlooked in the 
literature on 2D hole systems. Nevertheless, our main 
result in Eq. ([5]) bears some relation to earlier works. 
D'yakonov and Khaetskii2£ studied the Luttinger Hamil- 
tonian in an infinite square well and used the spherical 
approximation (72 = 73). They derived an expression 
for to || that agrees with our result in the appropriate 
limit. We also verified that Eqs. §5§ and ((BJ) can be ob- 
tained from a general k • p-approac h 26 i 27 after a lengthy 
calculation. In spite of the previous work, however, the 
relation of TO|| and g± to an additional parameter 7^ and 
the fact that 7^ is sensitive to F have been missing from 
the general knowledge of 2D hole systems. 

In conclusion, we showed that an external electric field 
can strongly modulate the perpendicular hole g-factor 
in SiGe QDs. By a detailed analysis, we ruled out the 
compositional-gradient mechanism as the origin of this 
electric-field effect. By analyzing the Luttinger Hamilto- 
nian in the 2D limit, we found a new correction term 7^ 
which had not been considered before in the literature. 
This new term, which corrects the "standard" expression 
for the HH g-factor, reflects the effect of a perpendicular 
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FIG. 5: Measured g-factor versus Vtg and Vb g for the trian- 
gular potential well configuration. The dispersion of the data 
points is too big to resolve a clear variation in the g-factor 
value. 

magnetic-field on the orbital motion, and it is ultimately 
related to the atomistic spin-orbit coupling of the valence 
band. 
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Appendix A: Supplementary data on the g-factor 
modulation 

According to the given theoretical description of the 
g-factor variation, we should be able to observe a strong 
effect just when the applied electric field drives the sys- 
tem close to the flat-band configuration. This is the case 
exhibited in Fig.|4]of the main text. In Fig.[5]thc g-factor 
variation is shown for the same device but in an electric 
field regime in which the potential well configuration is 
triangular. Although the range of Vt g and Vbg is similar 
in both Fig. 2] and Fig. [SJ there is no obvious change in 
the g-factor value for the latter case, as expected. 

It is known that thermal cycling can change the char- 
acteristics of a device and thus different set of data can 
be obtained for the same device^. A drastic change in 
Vtg (of about 3 V) appeared to produce a similar effect 
in our system. In this way, we could perform the g-factor 
measurements on two devices with completely different 
characteristics. This time, a slightly different type of 
measurement, which will be explained below, was per- 
formed in order to extract the g-factors. 

As discussed in the main text, by having a dual gate 



configuration it is possible to keep the same number 
of holes while changing the value of the external elec- 
tric field. The same Coulomb peak can be very easily 
followed since for each mV we move in Vt g , the peaks 
move ~ 13.2 mV in V& 3 . This ratio reflects the differ- 
ence in coupling between the QD and the two gates. 
Fig. [5] (a) shows how two coulomb peaks shift in Vb g 
while changing Vt g from approximately — 583 mV (blue 
trace) to V tg ~ — 620 mV (red trace). The reason why 
the Coulomb peaks appear to be split is that dI s d/dV s d 
was measured under V 3 d = 1 mV. 

By following particular Coulomb peaks we could, for 
different values of Vt g , sweep V s d as a function of Vbg, 
and obtain a stability diagram for the same hole state 
under different applied electric fields. From the Zccman- 
split lines observed in the stability diagrams, we have 
obtained the g-factors for 11 different values of Vt g and 
Vbg, see Fig. [5] (b). The inset which is surrounded by 
a blue rectangle shows the dI s d/dV s d as a function of 
(B,V s d) for Vtg — — 583.2 mV. This plot demonstrates 
that the parallel to the ground state line seen in the sta- 
bility diagrams insets (orange rectangles), is indeed the 
Zeeman splitting. From Fig. [5] (b) it becomes clear that 
the external electric field, has a very strong effect on the 
g-factor value. It reaches a minimum value of (2.3 ± 0.2) 
and a maximum value of (3.7 ± 0.4). This is, to our 
knowledge, the biggest g-factor value measured in SiGe 
nanostructurcs. Even larger g-factors can be obtained 
for QDs with higher Gc content, as has been shown by 
theoretical calculations^. 

Interestingly, two different behaviours are present in 
this set of data. From Vt g w —560 to Vt g w —630, the 
modulation in the g-factor resembles the one shown in 
the main text. It can be explained by the dependence of 
-fh on F. However, from V tg ~ —630 to V tg « —660 an 
abrupt increase of the g-factor for increasing electric field 
can be observed. This behaviour is not compatible with 
the explanation given in the main text. A careful inspec- 
tion of the stability diagram of Fig. H](b) for Vtg ~ —650 
shows that additional levels are present close the stud- 
ied level. No closely-lying levels are present in the range 
V tg G [—560,-630]. We believe that the data in the range 
V tg G [—630, —660] refers to a different size-quantization 
level than the data in the range V tg G [—560, —630]. In- 
deed, the ground-state level might have changed in going 
from one range to the other, which is presumably sig- 
nified by the pronounced cusp in the g-factor data at 
V tg f» — 630 mV in Fig. [HI (b). Such transitions are ex- 
pected to occur when the Fermi level is placed sufficiently 
deep in the valence band. Then, ladders of levels belong- 
ing to different heavy-hole subbands move with respect 
to each other when the electric field is varied. Since these 
levels have very different z-componcnts of the wave func- 
tion, they interact weakly and can come close to each 
other without a sizable level repulsion. Levels belonging 
to one ladder cannot come close to each other when the 
electric field is varied. We therefore speculate that the 
data in the range Vt g G [—630, —660] refers to a level 
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FIG. 6: (a) Plot of dI s d/dV s d versus Vbg for V a d = lmV show- 
ing how two coulomb peaks move while changing the value of 
the external electric field, (b) Plot of the g-factor vs Vbg and 
Vtg showing a non monotonic behavior of the g-factor value 
for different perpendicularly applied electric fields. The in- 
sets show stability diagrams at B — 1.5 T. Vbg is swept by 
50 mV while V s d by 2.4 mV. The blue inset demonstrates that 
indeed the parallel to the ground state line is due to the Zee- 
man splitting, because it merges with the ground state for 
vanishing magnetic fields. 



from the second heavy-hole subband. 

The following question arises: How can the two heavy- 
hole subbands have different characteristic values of the 
electric field at which the g-factor reaches maximum. It 
is important to note that our model in the main text is 
strongly simplified. The real confining potentials Vh(z) 
and Vi(z) in the NC are, most likely, never perfectly sym- 
metric and differ from each other. Therefore, it is natural 
to expect that the alignment along z of wave functions of 
the two heavy-hole subbands with respect to each other 
and, at the same time, with respect to the light-hole sub- 
bands is not perfect. The data in Fig. [5] (b) is consistent 
with the assumption that the wave function of the second 
heavy- hole subband is shifted towards the base of the NC 
and it begins to align with the light-hole wave functions 
at a later value of the electric field. 

It is important to remark that higher-order corrections 
to the expansion of the g-factor of the top-most subband 
cannot explain the abrupt increase taking place for Vt g € 
[—630, —660], because their contribution is small by ~ 
(w/d) 2 , which amounts to only a 10%-correction to the 
main term for our devices. 



Appendix B: The 2D limit of the Luttinger 
Hamiltonian 



Here we derive the expression for jh given in Eq. (6) 
of the main text. Our derivation is based on an ex- 
pansion of the Luttinger Hamiltonian around the two- 
dimensional (2D) limit, which we have recently outlined 
in R ef. [23l We summarize briefly the relevant results of 
Ref. [2J and address the following question: What is the 
effective Hamiltonian of the top-most hole subband in 
the 2D limit? 

We start with the representation of the Luttinger 
Hamiltonian in block form, 



H 



Huh Hm 
Hih Hu 



(Bl) 



where H^h and Hu are the main blocks, describing heavy 
holes and light holes, respectively. The off-diagonal 
blocks Hhi and Hi% = (Hhi) give the mixing between 
the heavy-hole and light-hole sectors and are responsible 
for avoided crossings which occur between heavy-hole and 
light-hole branches at higher energies. The energy axis is 
chosen to point downwards for holes. In the 2D limit, the 
heavy-hole and light-hole sectors become well separated 
in energy and the top-most hole subband (lowest in en- 
ergy) is described by purely heavy-hole states, i.e. states 
containing no admixture from the light-hole sector. It is, 
therefore, customary^ to neglect the off-diagonal blocks 
in Eq. (|BT[) when considering the 2D limit and describe 
the top-most hole subband by an effective Hamiltonian 
derived only from the block H^h, 



Ho 



(H 



111: I 



(B2) 



where (. . .) stands for averaging over the motion along z. 
Contrary to the common expectation, Eq. (|B2[) does not 
describe correctly the 2D limit of the Luttinger Hamil- 
tonian. The correct expression in the place of Eq. (|B2[) 
involves terms of the order of H^/Hu, which remain fi- 
nite in the 2D limit. While these terms do not change 
the qualitative picture of the result, they are important 
when mechanisms of electric control of the g-factor and 
effective mass are considered. 

Before proceeding further we need to clarify what ex- 
actly we mean by the 2D limit. By the 2D limit of the 
Luttinger Hamiltonian, we mean a limit in which the z- 
sizes of both heavy-hole and light-hole wave functions 
are sent to small values, e.g. by means of confining the 
holes to a thin layer of material. Of course, the Lut- 
tinger Hamiltonian is by itself an effective Hamiltonian 
(intended for top of valence band) and it is valid only for 
energies much smaller than the spin-orbital energy Aso- 
However, Aso is, typically, a large energy and one can 
envision the following limit, 



E < A z < A 



so, 



(B3) 



where E is the energy measured away from the edge of the 
top-most subband and A z ~ (72/^1) (k 2 ^ is the splitting 
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energy between heavy and light holes due to confinement 
(or due to strain if strain dominates). We, therefore, 
have in mind Eq. (|B3j) when talking about the 2D limit 
of the Luttiger Hamiltonian. It is important to note that, 
in the Luttinger Hamiltonian^, Aso is already sent to 
infinity and, thus, the second inequality in Eq. (|B3j) is, 
formally, fulfilled within the model. In practice, however, 
our derivation remains qualitatively correct up to A 2 < 
Aso- For Ge, one has Aso ~ 0.3 eV. 

After expanding the Luttinger Hamiltonian in terms 
of the small parameter w/d <C 1 (see Ref. [23h , the main 
blocks in Eq. (|B1[) are given at the leading order by 



H, 



hi, 



,2\ , 71 ~ 2 72 2 



71+72 ,,2 

+^bt -g h -B + U(x,y) + V h {z), 
7i - 72 / 



2m 
1 



k 2 ) I 7l+2 72 /; 2 

L y )+ 2m ** 



l i B (T-g l -B + U(x,y) + Vi(z), (B4) 



where k x and k y are the 2D momentum operators, which 
contain only the z-component of the magnetic field. The 
inplanc components of the magnetic field contribute to H 
at higher orders of the expansion, but not at the leading 
order considered here. The g- factors, gh and gi, entering 
Eq. (IB4[) , arc diagonal in the frame (x, y, z) and given by 



gh 







-6k- %q 



(B5) 



and 



91 



4k + 10q 
4k + Wq 
2k + 



(B6) 



The confining potential is assumed to separate into an 
inplanc component J7(a;, y) and a transverse Vh/i(z). The 
confining potential can be different for heavy and light 
holes, because of the strain. 

The off-diagonal blocks are related to each other by 
hermiticity, 



H M = (Hi h ) ] . 
For Hih, we keep the leading-order term 



Hi h 



{k X CTy + ky(T x ) k z 



(B7) 



(B8) 



In writing Eq. (|B4[) and (|B8[) . we made a choice of basis 
for the hole pseudo-spin. We chose the basis such that the 
pseudo-spin transforms under the time-reversal operation 
as a spin 1/2, see Ref. [H, 

|t>, = |3/2,-3/2), ||) fc = |3/2, +3/2), 

|t), = |3/2,+l/2), |4-) ; = 13/2,-1/2). (B9) 



This choice of basis is convenient since it allows us to 
have "standard" time-reversal transformations for all op- 
erators. In Eq. (|B9[) . we use basis of states of the angular 
momentum J = 3/2^1: 



= -T2 {x + tY)t 

= ±=(X-iY)i, 

= -j= [- (X + iY) 1 +2Z t] , 



V6 



[(X - iY) t +2Z |] , (BIO) 



where the functions X, Y, and Z are real and represent 
the Bloch amplitudes of the valence band in the absence 
of spin-orbit interaction. 

Next, we rotate away the off-diagonal blocks in 
Eq. (|B1|) using perturbation theory^ 2 -. We are allowed to 
do so if we focus on a heavy hole state which is far away 
from any light hole states. This requirement is fulfilled in 
an energy window between the edges of the first heavy- 
hole and light-hole subbands, i.e. close to the top of the 
valence band. Of course, occasionally, such a premise 
may hold for a state deep in the valence band. And, in 
particular, it may hold for a light-hole state. In the lat- 
ter case, we would be deriving an effective Hamiltonian 
in the favour of the light hole. However, the majority of 
states deep in the valence band do not permit application 
of perturbation theory. 

In order to quantify the applicability of our method, 
we introduce characteristic values for the momenta oper- 
ators in a localized state: 



ka = V( h \ k l (a = x,y,z), 



(Bll) 



where \h) stands for the localized (heavy-hole) state. 
Then, the expansion around the 2D limit (i.e. separa- 
tion of variables) is valid if 



(B12) 



whereas the perturbation theory can be applied if 



73 

m 



Uh Ll 



Eh, 



V/, 



(B13) 



where the index I refers to the light-hole states. Here, 
Eh and Ei are energies of the heavy-hole and light-hole 
states, respectively. Since we restrict our consideration 
to small inplanc momenta, we can approximate Eh and 
Ei by 



E h 



7i 



2 72 /)2N 



m 

7i + 272 



(k 



z/h> 



2 ) 

z/l' 



(B14) 



where (k 2 } h /i is the average of k 2 with the z-component 
of the wave function. With this approximation, the en- 
ergy denominator in expressions of perturbation theory 
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reads 



m v z 



^(kx&y H~ kyC 



(kl 



lh) ■ 
(B15) 

It is important to remark here that this energy denom- 
inator does not depend on the quantum numbers of the 
inplanc motion. To simplify the notations above, we de- 
noted all the quantum numbers by a single index, h (or 
I). Below, we shall single out the quantum number refer- 
ring to the z-componcnt of the wave function. Thus, we 
shall replace E t - E h in Eq. (fHTg]) by E l n - E% Q , where 
n and no denote quantum numbers of the motion along 
z and the superscripts h and / indicate to which type of 
hole the expression belongs. The z-components of wave 
functions associated with subband energies E l n and E% 
are denoted by ip l n (z) and tp^ (z), respectively. 
We average over the motion along z and obtain 



Next, note that 



— k„ 



i [k x , ky] 



Since [k x , k y 



-i(he/c)B z , we obtain 



lh,n 

2m 



(B17) 



(B18) 



(B19) 



where 



lh,n a — 



673 



E 



E l 



E h 

n 



(B20) 



Elctt = 



(H hh ) + (AH hh ), 



E l - E h 

^n no 



Therefore, we obtained a g-factor renormalization and 
(ip l J Hih \ip h ) , a mass renormalization. The heavy- hole g-factor g z 
changes to 

(B16) 

27 



where (...) = ■ • • | ^Pno) an d n o denotes the number 

of the heavy-hole subband for which we derive the effec- 
tive 2D Hamiltonian. In practice, hq can be one of the 
several first subbands and in the main text we set uq = I. 
Further, the derivation continues as follows, 



-6k 



The inplane mass reads 



m 



(AH hh ) = -^Ma 



71 + 72 - lh,: 



(B21) 



(B22) 



El 



E h 

^n a 
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